Shear viscosity and the r-mode instabihty window in superfluid neutron stars 



Cristina Manuel^'B and Laura 

^Instituto de Ciencias del Espacio (lEEC/CSIC) Campus Universitat Autdnoma de Barcelona, 
Facultat de Ciencies, Torre C5, E-08193 Bellaterra (Barcelona), Catalonia, Spain 
^Frankfurt Institute for Advanced Studies, Johann Wolfgang Goethe University, Ruth-Moufang-Str. 1, 60438 Frankfurt am Mam 

(Dated: December 11, 2012) 

We analyze how recent computations of the shear viscosity rj in the core of superfluid neutron stars 
affect the r-mode instability window. We take into account the contribution of superfluid phonons 
to the viscosity, both in their hydrodynamical and ballistic regime. Moreover, we also consider the 
recent computation of rj arising from the collisions of electrons with electrons and protons dominated 
by transverse plasmon exchange. We also discuss how the interactions among superfluid phonons 
and electrons might contribute to the shear viscosity, but argue that these interactions will be not 
be relevant for the study of r-mode damping. At very low temperatures, T < 10* K, the electron 
contribution to r; dominates, as then phonons are in a ballistic regime, and they mainly interact 
with the crust of the star. At higher temperatures the superfluid phonon contribution to 77 is purely 
hydrodynamic and starts to dominate the process of r-mode damping. While our results for the 
instability window are preliminary, as other dissipative processes should be taken into account as 
well, they differ substantially from previous evaluations of the r-mode damping due to the shear 
viscosity in superfluid neutron stars. 

PACS numbers: 04.40.Dg,97.60. Jd,26.60.-c,97.10.Sj 



I. INTRODUCTION 



Neutron stars are some of the most curious stellar objects in the Universe. They typically have a mass of approxi- 
mately 1AM sol and a radius of about 10 Km, and supranuclear densities are thus expected in their cores. Our present 
understanding of the nuclear interactions at finite density suggests that neutron stars are composed by well-defined 
diS'erent layers 0,01 • Thus, they are believed to have a crust of about 1 Km, with an outer part made of a lattice of 
ions embedded in a liquid of electrons and an inner region made of nuclei embedded in a liquid both of electrons and 
^5*0 superfiuid neutrons. In the interior of the star nuclei are melted and both neutrons and protons are expected to 
condense into BCS-like superfiuids. However, the neutron interaction in the ^S'o state at supranuclear matter density 
is repulsive, but it is still possible to form Cooper pairs in the '^P2 channel. The proton density is much smaller than 
the neutron density, and thus one expects proton pairing in the isotropic ^S'o channel, at least in some part of the 
core of the star. It is still unclear the content of the inner core of the star, as it could be formed either by hyperonic 
matter, phases with pion or kaon condensates, deconfined quark matter in a color superconducting phase, or other 
exotic dense matter phases [1]. 

One possible way to verify the internal structure of a neutron star is through the observation of the spectra of 
the different stellar oscillations, in a similar way as asteroseismology is used to reveal the internal structure of the 
Sun or other stars. In a neutron star there are different types of hydrodynamical oscillation modes. Of particular 
relevance are the r-modes, which are toroidal modes which only occur in rotating stars, the Coriolis force acting as 
their restoring force. R-modes are generically unstable in all rotating stars through their coupling to gravitational 
radiation (GR) emission When dissipative phenomena damp these r-modes the star can rotate without losing 

angular momentum to GR. If dissipative phenomena are not strong enough, these oscillations will grow exponentially 
and the star will keep slowing down until some dissipation or non-linear mechanism can damp the r-modes. As most 
of the known neutron stars are pulsars, and their rotation frequencies are known with a lot of accuracy, it has been 
recognized that the study of r-modes could effectively be used to constrain the internal stellar structure. 

R-mode oscillations have been studied extensively in the literature and various damping mechanisms have been 
proposed When viscous damping of the r-mode is taken into account the star is stable at low frequencies, or 
at very low or high temperatures, but there is typically an instability region at enough high frequencies The 
specific details of this window depend on the model considered for the star. Since the spin rate of several neutron 
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stars falls in this instability window, most of the studies of the r-modes are focused in detecting particular damping 
mechanisms that may modify or eliminate the instability window. In particular, it has been claimed that surface 
rubbing between the core and the crust of the star [1, |9[ results in a viscous boundary layer which damps r-mode 
oscillations for sufficiently small frequencies. Other damping mechanisms have been studied in the literature, as for 
example, the effect of mutual friction in superfluid neutron star (lol - [l2j . It has also been recognized that if the r-mode 
amplitude gets sufficiently large, its growth should be stopped by non- linear dynamical effects [l3l - [l7l |. 

It is our goal to review the form of the r-mode instability window for neutron stars with superfluid cores. Our main 
motivation comes from the fact that in previous studies of r-modcs the contribution to transport phenomena arising 
from superfluid phonons had been ignored, so it was impossible to assess its relevance for the r-mode damping. All 
previous studies of r-mode oscillations of superfluid neutron stars are based in the hydrodynamical equations of two 
fluid motions, one corresponding to the neutrons, the other one to the charge-neutral fluid, made up by protons locked 
to electrons. When finite temperature effects are taken into account one additional fluid motion has to be considered. 
It is well-known that superfluid systems like ^He have additional degrees of freedom, and that a two-fluid description 
of the superfluid is necessary to explain phenomena such as the second sound associated to thermal waves. At low 
T this second fluid gives account of the dynamics of some collective modes, the superfluid phonons. The superfluid 
phonon is a Goldstone mode which appears due to the fact that the neutron condensate, regardless whether ^Sq or 
neutron pairing is considered, spontaneously breaks the U(l) baryonic symmetry. The contribution of superfluid 
phonons to different transport coefficients has been computed for compact stars made of superfluid quark matter in 
the color-flavor locked (CFL) phase [l8l - l2]| . together with their effect on the damping of r-modes [22|. However, a 
similar analysis has not been yet carried out for superfluid neutron stars. 

In this article we consider a simplifled model of neutron star made up by neutrons, protons and electrons, using 
causal parametrization of the Akmal, Pandharipande and Ravcnhall [231 (APR for short) equation of state (EoS) 
to describe the /3-stable nuclear matter inside the star. This EoS is a common benchmark for a nucleonic equation 
of state, which is widely used in neutron star matter calculations. In the core of the star, neutrons pair flrst in a 
^5*0 channel, up to values of the particle density n ^ no = 0.16 fm~^, the normal saturation density [2J]. Up from 
this value, we consider that neutrons pair in the channel, all the way up to center of the star. Thus, in this 
work we don't consider the possibility of exotic forms of dense matter in the inner core of the star. We assume a 
critical temperature for the transition from superfluid to the normal phase in the range of Tc ~ 10^*^ K. On the 
other hand, protons pair in the ^Sq channel, up to critical densities of tic ^ (2.3 — 2.6) riQ 0. Thus, electromagnetic 
superconductivity on ly p revails up to that critical value of the particle density. Values of the ^5*0 proton gap are 



model dependent [25|-|33j. The result of Ref. [31| suggests that this pairing survives up to n ~ 2.3 no with a maximum 



value of ^ 1 MeV for n ~ 1.3 no . For deflniteness, we use the model of proton pairing of Ref. [3l[ in our computations. 
However, we don't need the precise numerical value of the neutron gap. 

In this paper we make a flrst approach to evaluate whether the viscous damping associated to superfluid phonons 
may affect the r-mode instability window of superfluid neutron stars, following our recent computation of Ref. S^] ■ We 
leave for the future the study of how other dissipative effects associated to the superfluid phonons, such as including 
their contribution to the bulk viscosities coefficients or their impact on mutual friction [35|, might also modify the 
r-modc instability window. In analyzing the r-mode instability window we also take into account the results for the 
electron contribution to the shear viscosity of Ref. [36j . where it has been claimed that previous results of rye had 
clearly overestimated its value. We also take into account the dissipation that might arise in the ballistic regime of 
phonons, and incorporate it in an effective shear viscosity coefficient, as done for other superfluids, such as "^He. With 
all the above considerations we see that the r-mode instability window as arising only from shear viscous damping 
is modifled substantially from previous estimates. For a comparison with previous results with the same sort of 
simplifled model for superfluid neutron stars, although with the use of a different EoS, see for example Ref. p^ . 

This paper is structured as follows. In Sec. |lT] we review different computations of the shear viscosity in the 
superfluid core of neutron stars, first as arising from electron collisions with electrons or protons (Sec. Ill A| . and then 
due to the phonon interactions with the boundary and among themselves (Sec. Ill Bp . While the computation of ?/ 
of Ref. (3^ was done under the assumption of neutron pairing in a ^S'o channel, we argue here that those results 
can as well be used for neutron pairing in a '^P2 channel, and explain how to incorporate the ballistic regime of the 
phonons in an effective shear viscosity coefficient. In Sec. Ill CI we discuss on how the clectron-superfluid phonon 
interactions might contribute to the shear viscosity, although we estimate that these contributions should not be the 
dominant ones for the r-mode damping. We analyze in Sec. IIIII how the values just presented in Sec. |TT] affect the 
r-mode instability window, modifying it with respect to previous results. We present a summary of our results and our 
conclusions in Sec. IIVI In Appendix]^ we show the computation of the coupling constant of the electron-supcrfluid 
phonon interaction in the low density limit. 

We use natural units 7i = c = /cb = 1 in all intermediate equations unless otherwise stated. However we present 
numerical values of the shear viscosity, temperatures and lengths in S.I. units. 
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II. SHEAR VISCOSITY IN THE CORE OF SUPERFLUID NEUTRON STARS 



Knowledge of transport coefficients in the core of neutron stars is necessary in order to determine the possible 
damping mechanisms of different oscillation modes of the star. In this Section we review the computations of the 
shear viscosity in the core of superfluid neutron stars, as due to electron collisions with electrons or protons (Sec. Ill Al) . 
and as due to phonons interacting among themselves or with the boundary, see Sec. IIIBI We also study the superfluid 
phonon-electron interactions in Sec. Ill C[ and how these interactions might contribute to the shear viscosity in the 
system, although we argue that they might be irrelevant for the study of the r-mode damping. The values of the shear 
viscosity presented in this Section will be used in Sec. Illll to assess on how they affect the r-mode instability window. 



A. Shear viscosity due to electron collisions 

Most of the transport coefficients needed for the study of the stellar evolution and structure of neutron stars have 
been originally studied by Flowers and Itoh in Ref. [s^l ■ These authors evaluated the transport coefficients in the core 
of neutron stars, both assuming normal and superfluid phases. In the last case they noted that the shear viscosity was 
mainly due to electron collisions with electrons or protons. The pure nucleon contributions were, in turn, drastically 
suppressed. In their computation Flowers and Itoh assumed that the collisions were dominated by low momentum 
transfers, and made the approximation of assuming static longitudinal plasmon exchange. According to these authors 
the value of the viscosity can be expressed in the extreme relativistic limit as [s^ 
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where Pe,Pp are the Fermi momentum of the electrons and protons, respectively, m is the proton rest mass, and a 
is the electromagnetic fine structure constant. Using the typical values of the Fermi momentum of the electrons and 
protons in the core of the neutron star. Cutler and Lindblom realized that this expression could be well approximated 
by the simple formula 

rie = 6.0 X lOV T-^—^ , (2) 
cm • s 

where p is the (local) mass density in g/cm^ and T the temperature in Kelvin. This is the formula that has been 
typically used to estimate the r-mode instability window in neutron stars with superfluid cores. 

Recently Shternin and Yakovlev [s^ have reviewed all the computations of the electron shear viscosity in neutron 
star cores, both in superfluid and normal phases. These authors realized that in the computations of transport 
coefficients performed by Flowers and Itoh the role of the charged particle collisions due to the exchange of transverse 
plasmons had been neglected. The transverse plasmon contribution enhances the coUisional rates with respect to 
the computation where only static longitudinal plasmons are considered, and this leads to a reduction of the shear 
viscosity. This fact had been first recognized by Heiselberg and Pethick, in the study of the transport coefficients for 
a degenerate quark plasma (39j . In the absence of neutron superfluidity and proton superconductivity, the electron 
contribution is the leading contribution to t], and according to the authors of Ref. [13], this can be expressed as 

r]e = 4.26 X 10-26(;^ „)14/9r-5/3_^ ^ (3) 

cm • s 

where n is the baryon particle density in cm^"^, Xp is the proton fraction and the temperature T in Kelvin. Recently, 
this corrected value of the shear viscosity has been used to analyzed the r-mode instability in neutron stars with non- 
superfluid, non-superconducting cores [40l . [4l| , as it was realized that also the electron contribution to rj dominates 
over the baryonic one in the normal phase of nuclear matter [3^ . 

Although the presence of neutron pairing does not affect directly the value of rje, proton pairing does. This is so as 
protons also contribute to the photon polarization tensors. Also in this case, Shternin and Yakovlev [s^ found out 
that r]e is dominated by the collisions mediated by transverse plasmons, which lead to the value 

where Ue is the electron particle density and A is the value of the ^5*0 gap for the protons. 

In this work, in order to obtain the energy density, electron density, proton fraction, and proton and electron Fermi 
momenta, we use the APR equation of state in a causal paramctrization with parameters S = 0.2 and 7 = 0.6, as seen 
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FIG. 1: (Color online) Shear viscosity due to electron collisions assuming prevalence of longitudinal plasmon exhange [Eq. ((2|] 
and electron scattering assuming prevalence of transverse plasmon exchange. In the last case, proton pairing is taken into 
account up to some critical value of the density [Eq.Q], and no proton pairing above that value [Eq. ((3|]. The plots are given 
as a function of the particle density in units of normal saturation density, no = 0.16 fm"'', for three different temperatures, 
T = 10*K (sohd lines), T = lO^K (dashed hues) and T = 5 x lO'^K (dotted lines). 



in Rcf. 3- With regard to the ^5*0 gap for protons, results are model dependent (25l - [33| . The works |29|-[3H suggest a 
maximum value of approximately 1 McV for densities of '--^ 1.3 no while it survives up to n '--^ 2.3 tiq |31| . although the 
inclusion of three-body forces, self-energy effects and the polarization interaction might reduce the domain of pairing 
[3^ . In this paper we use the estimate of Ref. [3l| . 

In Fig. [1] we show the contribution to the shear viscosity given by electron scattering with electrons and protons (in 
the following we will refer them solely as electron collisions) as a function of the particle density in units of normal 
saturation density, uq = 0.16 fm^"^, for three different temperatures, T = lO^K (solid lines), T = lO^K (dashed 
lines) and T = 5 x lO^K (dotted lines). A particle density of uq = 0.16 fm^"^ corresponds to a mass density of 
po ~ 2.8 X 10^* g/cnr^. A temperature of T ~ 5 x lO^K is close to the expected transition temperature to a normal 
fluid in case of proton superconductivity (Tc ~ 0.57 A). We present the predictions of rje by Flowers and Itoh, which 
assumed dominance of longitudinal plasmon exchange, and which are expressed in Eq. and also the predictions 
of Shtermin and Yakovlev, which assumed that rje is dominated by transverse plasmon exchange. In the last case, we 
assume that proton pairing only prevails for densities in the range n < ric 2.3 riQ. Thus, we use the values of Eq. (j4]) 
up to that critical value, and use Eq. ([3]) for higher densities, where no proton pairing occurs. Note that there is a 
clear discontinuity in ?/e in this case, but this is an artifact arising from the fact that Eq. @ is not valid in the regime 
where the proton gap vanishes, A — >■ 0. The two curves should interpolate continuously, however the formula for T^e 
in the case of almost vanishing A is not known. Because the discontinuity only affects the value of 77 for the critical 
value of the density, then it won't affect much our analysis of the impact of rj on the r-mode instability window, where 
one integrates the values of ry for all the densities in the star. 

We observe that, in concordance with the statements of Shternin and Yakovlev, the old values of the shear viscosity 
used in almost all studies of r-modes have been overstimated by a few orders of magnitude. We also observe that 
proton pairing leads to a decrease of the shear viscosity, as one could naturally expect. 



B. Shear viscosity due to phonons in superfluid neutrons stars including finite size effects 

The phonon contribution to the shear viscosity, assuming neutron pairing in a ^S'o channel, has been recently 
computed in Ref. [11]. We review here the results of that reference, and further argue that they can be extended for 
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phases with neutron pairing, where superfluid phonons also exist. We also comment how finite size corrections 
have to be taken into account when the phonon mean free path becomes larger than the size of the core of the star, 
and thus phonons are in a ballistic, rather than hydrodynamical, regime. 

The basic phonon interaction rates can be derived using effective field theory techniques. The effective phonon 
Lagrangian is presented as an expansion in derivatives of the phonon field, the terms of this expansion being restricted 
by symmetry considerations. The coefhcients of the Lagrangian can be computed from the microscopic theory, through 
a standard matching procedure, and thus they depend on the short range physics of the system under consideration. 
The leading order phonon Lagrangian in a derivative expansion turns out to be universal for all those superfluids 
sharing the same global symmetries. It reads [1^ 



CIS = P{X) , X = - dt^ - , (5) 



where and ^„ are the pressure and neutron chemical potential, respectively, of the superfluid at T = 0, and 

tf is the phonon field. We thus see that the leading phonon interactions depend on the equation of state of the star. 
The phonon contribution to the shear viscosity rj turns out to be dominated by binary collisions, and it was found 
that, as it happens in other superfluid systems, ri oc l/T^, the proportionality factor depending on the EoS of the 



star. For the particular computations of Rcf. [3J] the equation of state for /3-stable nuclear matter obtained by APR 
[23j . in a causal parametrization introduced by Heiselberg and Hjorth- Jensen 0, was used. The APR EoS ignores 
the effect of neutron pairing. However, because neutron pairing only affects those neutrons which are closed to their 
Fermi surface, one does not expect that the effect of pairing might have a big impact in the EoS. This effect should 
be of the order where A is the value of the s-wave gap. Given the fact that the neutron gap is believed to reach 



values A ~ 3 MeV (although correlations tend to reduce its value l2J]) while ^„ ^ 900 MeV, this seems to be a good 
approximation . 

We now argue that the results of Ref. for the contribution of the superfluid phonons to -q are as well valid for 
the anisotropic neutron pairing, if one also ignores corrections of order where A is the angular average of 

the gap function. At this order in the approximation the breaking of rotational invariance due to the existence 
of an anisotropic condensate is not felt by the superfluid phonons. This fact has been explicitly checked in Ref. (isj 
for ■^P2 neutron pairing, after computing the superfluid phonon kinetic term, integrating out the neutron degrees of 
freedom, and realizing that at that order one gets the same results that in a rotational invariant (s-wave) superfluid. 
However, one should be aware of a relevant difference between the s-wave and p-wave superfluid phases. In the last 
case there is also a spontaneous symmetry breaking of the rotational symmetry, and in this case there are additional 
massless collective excitations. The Goldstone modes associated to the breaking of rotational invariance, the so-called 
angulons, should as well contribute to the shear viscosity. We leave that analysis for a future project. In fact, the 
anisotropy in the hydrodynamical equations describing the p-wave superfluid phases of the core of neutron stars has 
not been considered so far in the study of any of its possible oscillation modes. Our results for rj when there is 
neutron pairing should be considered as a lower bound in this case. 

The phonon contribution to rj found out in Ref. [s^] diverges at sufficiently low temperatures. However, this is an 
unphysical result, which does not take into consideration that finite size effects actually prevent the viscosity from 
increasing indefinitely at low temperature. In Ref. [s^l it was actually shown that at sufficiently low T the phonon 
mean free path typically exceeds the value of the radius of the star. Let us define the viscous mean free path by the 
formula 

??bulk = ^PphCs^ph , (6) 

where jybuik is the phonon contribution to the shear viscosity in the hydrodynamical regime, = J^^^ is the phonon 
contribution to the mass density, and is the speed of sound. For hydrodynamics to be valid, one demands the 
Knudsen number Kn to be small. The Knudsen number is defined as the ratio of mean free path versus the typical 
macroscopic scale, in our case, the radius of the core of the star Rc, thus A'„ = Following Fig. 3 of Ref. [s^, the 
Knudsen number becomes much larger than 1 for values of T ^ 10* K (the precise threshold T depends strongly on 
the value of the nucleon particle density one considers) . For Kn > 1 an hydrodynamic description of the phonons is 
impossible, and one should write a Boltzmann equation to describe the phonon dynamics. 

When the phonon mean free path exceeds the size of the superfluid core of the star, the transport properties of 
the phonons are mainly governed by their interactions with the crust, which is the boundary of the superfluid region, 
rather than by their self-interactions. For simplicity, we will neglect the structure of the crust, and assume that the 
crust is not in a superfluid phase. This is a sort of approximation that has been typically done in all studies of the 
r-modes oscillations, as even it is considered that r-mode damping might be dominated by the rubbing effect of the 
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fluid with either a rigid or clastic crust [8|, |9j . We make here the same assumption in considering that the superfluid 
phonon motion is confined to the core of the star. A much more rigorous treatment should consider that the superfluid 
phonons of the core can penetrate into the inner crust, where neutrons also pair in a s-wave condensate. This would 
also imply to consider that a transfer of neutrons between the inner crust and outer core is indeed possible, which 
should have consequences in the r-mode evolution, through the so-called rocket effect [13] . In this article we will take 
the simpler assumption by considering the crust as a rigid boundary. If the superfluid phonon motion is confined to 
the core of the neutron star, then phonons reaching the crust can be either absorbed or they can be scattered back. If 
the phonons arc diffused, they will exert a shear stress on the boundary, in the same way as the particles in a rarefied 
gas produce shear stress when diffused by the walls of their container. 

In the case described above, when > 1, phonons are in a ballistic regime, where no hydrodynamical description 
of their behavior is possible. However, it has been observed experimentally that for ^He at very low temperature 
ballistic phonons can still efficiently damp the movement of immersed objects, see Refs. (45l - l47j . This damping can 
still be described by an effective or ballistic shear viscosity coefhcient defined by 



where d is the typical size of the oscillating object, and x is a coefficient that measures the probability that the phonons 
are diffused at the boundary with the object. This expression leads to excellent agreement with the experimental 
data for large Knudsen numbers [47| . 

In the intermediate regime between large or small Knudsen numbers, that is, between the ballistic and liydrodynamic 
regime, one can effectively calculate the shear viscosity according to 



This formula can be justified from the fact that the shear viscosity is proportional to a coUisional relaxation time, and 
that if one considers different processes occurring independently with different relaxation times, the total relaxation 
frequency is the sum of the partial relaxation frequencies. The same sort of formula has been used recently to describe 
the low T behavior of the experimental values of the shear viscosity in the superfiuid phase of a cold Fermi ato mic g as 
in the unitarity limit, where one assumes that it is also dominated by the dynamics of the superfluid phonons [48|. 



We will use Eq. (|8]) to give account of the phonon contribution to r] in the core of the star, taking for d the value 
of the radius of the core. The radius of the core is fixed by the transition between the crust and the core, which 
corresponds to n ^ 0.5 np. The value of x could be determined if at a microscopic level one could determine the sort 
of scattering processes that the phonons experience when colliding with the crust. Here we will assume that x ~ 1. 
We take for Tybuik the values of the shear viscosity computed in Ref. [s^l, which are strictly valid in the hydrodynamic, 
or low Knudsen number, regime. Note that Tyeff gets dominated by the smaller contribution among r^buik and Tybaii, so 
it describes properly the different hydrodynamical, transition and ballistic regimes. 

In Fig. [5] we show values of the ballistic shear viscosity 77baii (dotted lines), of the hydrodynamic shear viscosity 
Vhuik (dashed lines) and the effective viscosity ryeff (solid lines) that allows to describe all the regimes (ballistic, 
hydrodynamical, and the transition among the two). These contributions arc displayed as a function of temperature 
for different particle densities in units of the normal saturation density. The densities studied are those expected to 
be reached in the core of neutron stars. The ballistic contribution in Fig. [5] is calculated for a core radius Rc = 9.96 
Km, which corresponds to a density of 0.5 ng for a star of 1.93 Mgoi- 

We observe three distinct regions depending on the dominant contribution to the shear viscosity. For T < lO^K 
the shear viscosity is described in terms of the ballistic contribution for all densities in the core. As the temperature 
increases, the transition from the ballistic to the hydrodynamic domain takes place. This intermediate regime is 
strongly dependent on the density, as we can see in Fig. [21 For n ^ uq this transition domain appears for T ~ 
lO^K, while for higher temperatures, the hydrodynamic reg ime dominates. As indicated before, the validity of the 
hydrodynamical description was already studied in Ref. [3J] by means of the analysis of the phonon mean free path 
as compared to the radius of the neutron star. As one would expect, this study and that of Ref. [s^l are in good 
agreement. We also observe that for a fixed value of the temperature the ballistic viscosity decreases as soon as the 
particle density increases, while the hydrodynamical viscosity exhibits the opposite behavior. 

The question that remains is whether, in superfluid neutron stars, the phonon effective shear viscosity could govern 
the r-mode damping inside the core. As discussed in Sec. Ill Al it is believed that the electron contribution to the 
shear viscosity dominates with respect to any other baryonic contribution. In Fig. [3J we present the shear viscosity 
given by electron collisions assuming prevalence of transverse plasmon exchange [Eqs. ([3]) and Q] together with the 
effective phonon shear viscosity, where the ballistic and hydrodynamical regimes are both considered [Eq. These 
are shown as a function of the particle density for three different temperatures, T = lO^K (left panel), T = IQ^K 
(middle panel) and T = 5 x lO^K (right panel). 
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FIG. 2: (Color online) Phonon contributions to the shear viscosity as a function of temperature for different particle densities 
in units of normal saturation density, no = 0.16 fm^"^. We plot separately the ballistic viscosity jybaii (dotted lines), the 
hydrodynamical viscosity r^buik (dashed lines) and the effective description rjcs (solid lines), which takes properly into account 
which regime dominates the dissipation. 



We concentrate now on the analysis of the phonon shear viscosity and the comparison with the viscosity coming 
from electron collisions. We observe that the effective description of the phonon shear viscosity strongly depends on 
the temperature. For T ~ lO^K (left panel of Fig. [3]), the ballistic description dominates and the viscosity scales 
as 1/Cj with density. However, as the temperature increases, the transition between the ballistic and hydrodynamic 
domains takes place. For T ^ lO^K (middle panel of Fig. [3]), we find that the phonon shear viscosity shows a strong 
density dependence coming from this transition. As seen in Fig. [2J low densities reach faster the hydrodynamic 
domain. Thus, for T ~ lO^K and n ^ uq, the shear viscosity increases rapidly as function of the density following 
a hydrodynamical behavior. This is followed by a decrease with density as the ballistic description takes over. It is 
interesting to note that for T ~ lO^K and n ^ 2.5no, phonons interacting among themselves and with the boundary 
provide a higher dissipation than electron collisions. For T ~ 5 x lO^K (right panel of Fig. [3]), the description of 
the phonon shear viscosity is hydrodynamic and phonon-phonon scattering dominates over electron collisions in the 
whole core. Therefore, we conclude that processes involving phonons can have a dominant role in the dissipation in 
neutron stars for some values of the temperature, and affect, in turn, certain oscillation modes of the star, such as 
the r- modes. We will address this issue in Section Hill 



C. Interactions among superfluid phonons and electrons in the core of neutron stars 

We study here the microscopic interactions of the superfluid phonons with the other light degrees of freedom of 
the core of the neutron star, the electrons. This may allow us to see whether these processes can give a substantial 
contribution to the shear viscosity. 

Let us assume that both neutrons and protons form Cooper pairs in a channel (a similar reasoning can be applied 
for the p-wave superfluid case). Effective field theory techniques can be used to write down the Lagrangian describing 
the electron and phonon interactions. At leading order in a derivative expansion, simple symmetry arguments suggest 
that this is given by 

-Cph-cl = 9 (V'e7°V'e5o'^ + asi'el^'ipediip) , (9) 

where ^pe is the electron field. The values of the coupling constants g and can be obtained after integrating out 
the neutron, the proton and the phonon from the complete underlying microscopic field theory. In particular, their 
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FIG. 3: (Color online) Shear viscosity resulting from electron collisions assuming prevalence of transverse plasmon exchange 
[Eqs. ((3| and Q] together with the effective description of the phonon shear viscosity [Eq. ((8|], as a function of the particle 
density in units of normal saturation density, no = 0.16 fm^'*, for three different temperatures, T = 10*K (left panel), T = lO^K 
(middle panel) and T = 5 x lO^K (right panel). 
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FIG. 4: Feynman diagram for electron-phonon scattering in matter with protons and neutrons forming ^Sq Cooper pairs. 



values can be obtained in the low density domain after the evaluation of the Feynman diagram depicted in Fig. |4] (see 
Appendix lAl for details). At low density one finds 



2p„mcs 



(10) 



where p„ if the neutron Fermi momentum, p„ is the neutron mass density, and a„p is the neutron-proton scattering 
length, that in vacuum takes the value a„p = —23.7 fm. One might expect, however, that at high density the above 
value of g might be substantially corrected. We will assume that Os ^ 1. 

In the supcrfluid core of a neutron star the electron chemical potential has values of a few hundreds of MeV, 
while the temperature is very low, T <C 1 MeV. So it is a good approximation to consider that electrons are totally 
degenerate, and neglect how thermal effects may modify their Fermi sea. The interactions of the phonons with the 
electrons inside the Fermi sea are very suppressed, and a phonon. which typically carries an energy of order T with 
T <^ pLf., can hardly excite an electron inside the Fermi sea. However, we note that phonons can interact with the 
electrons which lie on their Fermi surface. 
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Let us see now how knowledge of Quantum Electrodynamics (QED) in the ultradegeneratc limit can be used to infer 
properties of the superfluid phonon-electron interactions. For that purpose, it is enough to realize that the lagrangian 
Eq. ([9]) can also be written as 

'Cph-ol = g^'el^lpeA^ , (11) 

if we define = (dgip, agdiip). In this way, this interaction term resembles that of the gauge interaction of the 
electrons with an electromagnetic field potential. The phonon self-energy can thus be obtained from the knowledge 
of the photon self-energy. In momentum space, with wave vector = (fco,k), one has 

\A^{-k)^^^^^{k)K{k) - -\p{~k)k^^^^^^{k)K^,(k) , (12) 

where we have defined the wave vector A;^ = (fco, a^k), and we take the QED expressions for the polarization tensor 
simply replacing the electromagnetic coupling constant e by 5. Expressions for the polarization of QED in the limit 
of high chemical potential /^e and vanishing temperature are well-known for values of the momenta much less than 
/Ltg. These are given by the so called hard dense loops (HDL), and it can be shown that only the electrons on their 
Fermi surface contribute to the polarization tensor [49|. We thus see that the superfluid phonon polarization tensor 
is expressed in terms of the photon polarization tensor according to 

n(fc) - -kjz,\\^^^^(k) . (13) 

The HDL polarization tensor contains an imaginary part which gives account of the so-called Landau damping 
[50| . These are processes where a photon is absorbed (or emitted) by the electrons. This a coUisionless dissipation 
which does not involve an increase of entropy, and thus such a decay does not contribute directly in the evaluation 
of the transport coefficients of the system. Based on the parallelism we have just traced of the photon and phonon 
interactions with the electrons, we are then lead to conclude that the phonons in the superfluid nuclear media also 
suffer Landau damping. The phonon damping rate, which is directly related to Landau damping, is given by 

^^ ImH(fco=i^.k) ^^,^,^^^^^^,^^^ (14) 

where we have used that the phonon dispersion law is Ek = Cgk at leading order. 

Because Landau damping does not involve an increase of entropy, it cannot contribute directly to the evaluation 
of the shear viscosity. Following the analogy with the QED case, one realizes that electron collisions mediated by 
one-phonon exchange will turn out to be the most relevant process involving both electrons and phonons that might 
contribute to the shear viscosity. We leave that computation for a future project, but present below some arguments 
of why our conclusions for the evolution of r-modes might not be affected by this process. The computation of the 
shear viscosity as due to the collisions mediated by one-photon exchange or one-phonon exchange in the degenerate 
plasma turn out to be very similar. For the QED case the scattering matrix is given by [39j 

M ~ e\j;D^^''{k)Jt , (15) 

where and are the electromagnetic currents associated to the two electrons participating in the collision, and 
D^^^ is the photon propagator. In the electron-phonon case the scattering matrix reads 

M ^ g'j;,k^^k''G{k)rJ , (16) 

where G{k) is the phonon propagator. Because the energy transfer in these collisions is at most of order T [s^, one 
is lead to conclude that the relevance of one process versus the other one is provided by the ratio e^/(.g^T^). We have 
evaluated this ratio, using for this purpose the value of g given in Eq. pUj) which is at most 0.1 MeV~^. We conclude 
that at temperatures T ^ (10^ — 10^)K, the electron collisions mediated by one-photon exchange are more relevant 
than those mediated by one-phonon exchange. 

For the reasons mentioned above, we will neglect in this paper the shear viscosity arising from the interactions 
among electrons and the superfluid phonons. A much more detailed discussion about the superfluid phonon-electron 
interactions will be presented somewhere else. 
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FIG. 5: (Color online) R-mode instability window for superfluid neutron stars. The critical (normalized) frequency is plotted 
for two neutron star mass configurations (1.4 Msoi and 1.93 Msoi) as a function of temperature for the different dissipative 
processes studied: electron collisions assuming prevalence of longitudinal plasmon exhange [Eq. ^] (dotted lines), electron 
scattering assuming prevalence of transverse plasmon exchange [Eqs. ((S]) and (|4])] (dashed lines), the effective description of the 
phonon shear viscosity (dashed-dotted lines), and the sum of both the shear viscosity due to electron collisions with prevalence 
of plasmon exchange and the effective phonon shear viscosity (solid lines). Here flo = y/Gnp, where p is the average mass 
density of the star. 



III. SHEAR VISCOSITY AND R-MODE INSTABILITY WINDOW 



In order to assess the r-mode instability window one has to compute the value of the different time scales associated 
to both the instability due to gravitational wave emission tgr , and to the various dissipative processes that may damp 
the r-mode in the star. In a superfluid neutron star the critical rotation rate is obtained by studying the equation 

1 _ 1 1 1 1 

where t,;,t^ and tmf refer to the damping times associated to processes that involve the shear viscosity, the bulk 
viscosity and mutual friction, respectively. For superfluid neutron stars the bulk viscosity damping time has been 
assumed to be negligible in the literature. Only the effects of the shear viscosity due to electron collisions, and the 
mutual friction processes arising from the scattering of the electrons off the magnetic fields entrapped in the core of 
the superfiuid neutron vortices have been taken into account (lol - [l^ . In this article we only analyze the effect of the 
shear viscosity in the r-mode instability window. The superfluid phonons also contribute both to the bulk viscosities, 
and to mutual friction, but the microscopic computations necessary to evaluate their effect on the damping of the 
r-modes have not been yet carried out. 

In this paper we simply study the impact of the shear viscosity in the r-mode instability window. We consider both 
the shear viscosity as arising from electrons and phonons, and take the values discussed in Sec. HI] [HI]. Thus we study 
the roots of the curve 

- + 77n = ' (18) 

where the characteristic time scales are given by the expressions of Ref. (in S.I. units) 



21+2 ^B. 

.21+2 



pr^'+^dr , (19) 
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and 



{l-l){2l + l) 7jr^'dr\ pr^'+^dr] . (20) 



We study only r- modes with I ~ 2 because they are the dominant ones 0, IB] ■ 

In Fig. [5] we show the r-mode instabihty window of the neutron star as arising from considering only shear viscous 
damping for two star mass configurations, 1.4 Mgoi and 1.93 Mgoi. The corresponding radius of those configurations 
are 11.7 Km and 10.2 Km, respectively, while the central densities are, in turn, 3.2 np and 7.6 ng. Let us recall again 
that we are using the APR equation of state in a causal parametrization of Ref. [21 in order to assess the value of 
the mass density as a function of r, as well as the values of the shear viscosity. The critical frequency in Fig. [S] is 
given in terms of Oq = sJOixp, where p is the average mass density of the star, and it is displayed as a function of the 
temperature for the different processes studied. We plot the old estimate assuming dominance of electron collisions 
with prevalence of longitudinal plasmon exhange [Eq. ^\ (dotted lines). We also plot the estimate that arises by 
considering electron scattering assuming prevalence of transverse plasmon exchange [Eqs. (O and (jH)] (dashed lines) 
and the effective description of the phonon shear viscosity (dashed-dotted lines). Moreover, we display the sum of 
both the shear viscosity due to electron collisions with prevalence of transverse plasmon exchange and the effective 
phonon shear viscosity (solid lines). 

We observe that for low temperatures T < 5 x lO^K for a star of 1.4 Mgoi, and T < lO^K for a mass configuration 
of 1.93 Msoi the electron contribution to the shear viscosity dominates and governs the processes of r-mode damping. 
However, the early computations of rje had clearly overestimated the value of r]e by some few orders of magnitudes, as 
shown in Fig. [U which lead also to an underestimation of the r-mode instability window. At very low T the critical 
frequency turns out to be almost two times smaller than that predicted by the old computation of rj by Flowers and 
Itoh. 

The dissipation stemming from phonon processes takes over for T ~ 5 x lO^K for a star of 1.4 Mgoi while temperatures 
of T ^ lO^K are needed for a mass configuration of 1.93 Mgoi- This can be understood as follows. A star of 1.4 Mgoi has 
a central density of n '--^ 3 ng. And, from the middle panel of Fig. |3l we see that for T ^ lO^K the processes involving 
phonons govern the dissipation for densities up to 2.5 no. Thus, for a low-mass star of 1.4 Msoi and T ^ lO^K, the 
phonon shear viscosity will be the dominant dissipative process. Indeed, as seen in Fig.[5l even for lower temperatures 
of T ^ 5 X lO^K this would be the case. For a mass configuration of 1.93 Mgoi, higher densities are reached at the 
center of the star of rt ~ (7-8) no . Then, higher temperatures are needed so that the processes involving phonons 
will govern the dissipation (middle and right panel of Fig. In this case, the temperature regime where phonon 
processes dominate over electron collisions is of T ^ lO^K . These results can be inferred directly from the solid lines 
in Fig. O From the above considerations we can conclude that phonon processes should be taken into account for 
T > (lO'^ — 10^)K in order to assess the r-mode instability window. 



IV. CONCLUSIONS 



In this article we have made a first attempt to assess whether the dissipation processes arising from the scattering of 
supcrfluid phonons can have any impact on the r-mode instability window associated to neutron stars with superfluid 
and superconducting cores. We have so far only considered how shear viscous damping in the core of the star affects 
the r-mode instability window, leaving for a future analysis a much more complete study. Up to now the r-mode 
instability window of superfluid neutron stars had been analyzed taking into account only the shear viscosity due to 
electron collisions and the mutual friction processes arising from the scattering of the electrons off the magnetic fields 
entrapped in the core of the superfluid neutron vortices, or the viscous Eckman layer effect. The superfluid phonons 
should also contribute to mutual friction, through their scattering with rotational vortices of the star (a computation 
of this effect in the ballistic regime of phonons was carried out for the CFL phase in Ref. [2l|), or to the bulk viscocity 
coefficients. But the microscopic computations necessary for the evaluation of their impact on the r-mode damping 
have not yet been carried out. 

We have considered a simplified model of a neutron star made up of neutrons, protons and electrons, described 
by a causal parametrization of the APR equation of state. The star is composed by a rigid crust, and a superfluid 
and superconducting core. While there are solid theoretical arguments to believe that neutron and proton pairing 
occur in the core of the star, the precise values of the different gaps are very model dependent. We assume a critical 
temperature for the onset of superfluidity in the range of Tc ~ 10^'' K. Our computations for the shear viscosity as 
arising from superfluid phonons do not depend on the value of the neutron gap, and then, we don't need to rely on 
any speciflc model for neutron pairing. We simply assume that superfluidity occurs in the whole core of the star. We 
consider one speciflc model for proton pairing, as the electron viscosity does depend on the values of the proton gap. 
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following the results of Ref. [3^ . In this way, the values for the electron viscosity presented in Fig. [T] are thus model 
dependent. However, given the similar values for the proton gap in (29l - [3l| . we expect that our results for the shear 
viscosity and the r-mode instability window will be not much affected by the model of proton pairing. 

We have considered two neutron star mass configurations, of 1.4 Mgoi and 1.93 Mgoi- The radius and central densities 
of these two sort of stars are quite different. Our analysis suggests that at very low T the electron shear viscosity 
dominates over the phonon viscosity when obtaining the r-mode instability window. Previous studies of the r-mode 
instability window were based on a computation of rje that was overestimating its value by few orders of magnitude, 
which results in an underestimation of the r-mode instability window. The recent computations of 7]^ of Ref. [36j 
suggests that critical values of the frequency can be up to half of the values that were predicted in the literature for 
sufficiently low T. The role of superfluid phonons starts to be relevant only at T ~ 5 x 10^ K for the 1.4 Mgoi star, and 
for T ~ 10^ K for the 1.93 Mgoi. This difference can be understood from the fact that the central densities reached in 
the two cases differ by almost a factor of two. At lower T superfluid phonons in the core are basically in a ballistic 
regime, and the dissipation associated to their interactions with the crust are not strong enough so as to compete with 
the dissipation associated with electron collisions. We thus find that superfluid phonons only contribute at relatively 
high T, and their effect might only be needed to study the r-mode evolution of hot young neutron stars. 

Since wc find a serious modification of the r-mode instability window for superfluid neutron stars as arising only 
from the shear viscous damping, the Eckman layer damping time scale for the r-modes should be also modiflcd. This 
is so because this effect, which measures the rubbing effect between the core and the crust of the star, also depends 
on the value of the shear viscosity in the core of the star. This analysis, together with a more complete study of all 
the dissipation mechanisms associated to superfluid phonons, will be the subject of future studies. 
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Appendix A: Evaluation of the coupling constant of the phonon-electron interactions 

The value of the coupling constants g and as of Eq. ([9]) can be computed from the microscopic theory if one integrates 
out the heavy fields of the theory. Assuming that the nuclcons are paired in a s-wave channel, the computation can 
be easily done in the low density domain. 

Let us start by recalling the microscopic theory that describes the nucleon dynamics when the momenta of the 
nucleons is small, i.e., smaller than the pion mass. The basic Lagrangian describing the nucleon interactions at low 
momentum can be described by the effective field theory [5l| 

(jPP fynp 

- -f- W(x)^p(x)) {i^lix)M^)) - -f- {,pUx)Mx)) {^Ux)M^)) +■■■ , (Al) 

where ipm'>Pp are the neutron and proton wave-functions, is the electromagnetic gauge potential, and Cq are 
coupling constants that can be related to the different nucleon-nucleon scattering lengths Uij, when these are small, 
Cq = 47ra.y/m. Note that wc have taken into account isospin breaking terms in the contact interactions, and that 
we have not displayed other contact interactions in a higher spin channel. 
It is convenient to write down the neutron wavefunction as 

i^n = ^n,oe"^ , (A2) 

where the phase of the neutron wavefunction is the Goldstone field. Then it is easy to check that the kinetic term 
associated to ipn^o is the same as that for the whole neutron wavefunction, where one only has to substitute the 
standard derivatives by covariant derivatives 

i^iiidt - ^ + ^^n)'^Pn = ipi^oi^idt + iAo) - ^^^^ M„)-0n,o , (A3) 
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where = {dt(p, 'V(p). Integrating out the field ipn,o, we obtain the Lagrangian for the superfluid phonon, see Eq. ([5|). 
The value of the coupling constants g and Og can be obtained after computing the diagram displayed in Fig. For 
the purpose of the present discussion, it will be enough to present a quick estimate of the coupling constant g. We 
do not calculate the value of as, as it turns out that it may not affect in any relevant way our conclusions about the 
phonon-electron interactions in the core of the star. The superfluid phonon momentum in the star is of the order of 
T/cs, where Cs is the value of the speed of sound, and it is much less than the Fermi momentum of the neutrons or 
protons in the star. Thus, one can evaluate the neutron and proton loops of the diagram of Fig. |3]by neglecting the 
phonon momentum in the loop. Every nucleon loop contributes as (m pp^N)/Tr'^, where pf.n is the Fermi momentum 
of the nucleon N (see for example Ref. (H^)- The photon propagator is also assumed to be dominated by the value of 
the Debye or Meissner photon mass. Thus, one concludes that 

2pF.nmCs ,,,, 
g ~ , — anp ■ (A4j 
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